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Abstract 

We give a concise axiomatic introduction into the fundamental structure of classical elec- 
trodynamics: It is based on electric charge conservation, the Lorentz force, magnetic flux 
conservation, and the existence of local and linear constitutive relations. The inhomogeneous 
Maxwell equations, expressed in terms of D l and Hi , turn out to be a consequence of electric 
charge conservation, whereas the homogeneous Maxwell equations, expressed in terms of Ei 
and B\ are derived from magnetic flux conservation and special relativity theory. The exci- 
tations D l and Hi , by means of constitutive relations, are linked to the field strengths Ei and 
B\ Eventually, we point out how this axiomatic approach is related to the framework of gauge 
field theory. 
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1 Introduction 



In nature one has, up to now, identified four fundamental interactions: Gravity, electromag- 
netism, weak interaction, and strong interaction. Gravity and electromagnetism manifest 
themselves on a macroscopic level. The weak and the strong interactions are generically 
microscopic in nature and require a quantum field theoretical description right from the begin- 
ning. 

The four interactions can be modeled individually. Thereby it is recognized that electro- 
magnetism has the simplest structure amongst these interactions. This simplicity is reflected 
in the Maxwell equations. They, together with a few additional assumptions, explain the elec- 
tromagnetic phenomena that we observe in nature or in laboratories. 

Without digressing to philosophy, one may wonder about the origin of the Maxwell equa- 
tions. Should we believe in them as such and just study their consequences? Or should 
we rather derive them from some deeper lying structures? Certainly, there are already some 
answers known to the last question. The Maxwell equations rely on conservation laws and 
symmetry principles that are also known from elementary particle physics, see ID ED- hi 
the framework of classical physics, authoritative accounts of electrodynamics are provided by 
|[T9ll23l . e.g.. In this paper we would like to add some new insight into this subject. 

We will provide a short layout of an axiomatic approach that allows to identify the basic 
ingredients that are necessary for formulating classical electrodynamics, see We believe 
that this axiomatic approach is not only characterized by simplicity and beauty, but is also of 
appreciable pedagogical value. The more clearly a structure is presented, the easier it is to 
memorize. Moreover, an understanding of how the fundamental electromagnetic quantities 
D\ Hi, E i: B % are related to each other may facilitate the formulation and solution of actual 
electromagnetic problems. 

As it is appropriate for an axiomatic approach, we will start from as few prerequisites as 
possible. What we will need is some elementary mathematical background that comprises 
differentiation and integration in the framework of tensor analysis in three-dimensional space. 
In particular, the concept of integration is necessary for introducing electromagnetic objects as 
integrands in a natural way. To this end, we will use a tensor notation in which the components 
of mathematical quantities are explicitly indicated by means of upper (contravariant) or lower 
(covariant) indices ll22l . The advantage of this notation is that it allows to represent geometric 
properties clearly. In this way, the electromagnetic objects become more transparent and can 
be discussed more easily. For the formalism of differential forms, which we recommend and 
which provides similar conceptual advantages, we refer to lfT0l l4l. 

We have compiled some mathematical material in the Appendix. Those who don't feel 
comfortable with some of the notation, may first want to have a look into the Appendix. Let 
us introduce the following conventions: 

• Partial derivatives with respect to a spatial coordinate x l (with i,j,--- = 1, 2, 3) or with 
resepct to time t are abbreviated according to 



_d_ 



di 



i -i 



d_ 

Of 



(1) 
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We use the "summation convention". It states that a summation sign can be omitted 
if the same index occurs both in a lower and an upper position. That is, we have, for 
example, the correspondence 



atP 1 . (2) 



i=i 

We define the Levi-Civita symbols and e ,jfc . They are antisymmetric with respect 
to all of their indices. Therefore, they vanish if two of their indices are equal. Their 
remaining components assume the values +1 or —1, depending on whether ijk is an 
even or an odd permutation of 123: 



1, for ijk = 123, 312, 231, 
-1, for ijk = 2\3, 321, 132. 



_ ijk _ ) ■*• > "J'" ' ~> ^ ' , 



With these conventions we obtain for the gradient of a function / the expression dif. The curJ 
of a (covariant) vector Vi is written according to e l ^ k djV^ and the divergence of a (contravari- 
ant) vector (density) w l is given by diW 1 . 

Now we are prepared to move on to the Maxwell theory. 



2 Essential classical electrodynamics based on four axioms 

In the next four subsections, we will base classical electrodynamics on electric charge con- 
servation (axiom 1), the Lorentz force (axiom 2), magnetic flux conservation (axiom 3), and 
the existence of constitutive relations (axiom 4). This represents the core of classical electro- 
dynamics: It results in the Maxwell equations together with the constitutive relations and the 
Lorentz force law. 

In order to complete electrodynamics, one can require two more axioms, which we only 
mention shortly (see J4| for a detailed discussion). One can specify the energy-momentum 
distribution of the electromagentic field (axiom 5) by means of its so-called energy-momentum 
tensor. This tensor yields the energy density (D l Ei + HiB l )/2 and the energy flux density 
e^ h EjH k (the Poynting vector), inter alia. Moreover, if one treats electromagnetic problems 
of materials in macrophysics, one needs a further axiom by means of which the total electric 
charge (and the current) is split (axiom 6) in a bound or material charge (and current), which 
is also conserved, and in a free or external charge (and current). This completes classical 
electrodynamics. 

2.1 Electric charge conservation (axiom 1) and the inhomogeneous 
Maxwell equations 

In classical electrodynamics, the electric charge is characterized by its density p. From a 
geometric point of view, the charge density p constitutes an integrand of a volume integral. 
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This geometric identification is natural since, by definition, integration of p over a three- 
dimensional volume V yields the total charge Q enclosed in this volume 



Q:= / pdv. (4) 



We note that, in the Si-system, electric charge is measured in units of "ampere times second" 
or coulomb, [Q] = As = C. Therefore the Si-unit of charge density pis [p] = As/m 3 = C/m 3 . 

It is instructive to invoke at this point the Poincare lemma. There are different explicit 
versions of this lemma. We use the form d69t that is displayed in Appendix lAl Then (if space 
fulfills suitable topological conditions) we can write the charge density p as the divergence of 
an integrand D l of a surface integral. Thus, 



dtD^p (divV = p). (5) 



This result already constitutes one inhomogeneous Maxwell equation, the Coulomb-Gauss 
law. In parenthesis we put the symbolic form of this equation. 

Electric charges often move. We represent this motion by a material velocity field u\ that 
is, we assign locally a velocity to each portion of charge in space. The product of electric 
charge density p and material velocity u 1 defines the electric current density J\ 

f = pu l . (6) 

Geometrically, the electric current density constitutes an integrand of surface integrals since 
integration of J 1 over a two-dimensional surface S yields the total electric current / that 
crosses this surface, 

f dai . (7) 

5 

We have, in Si-units, [/] = A and [J*] = A/m 2 . 

We now turn to electric charge conservation, the first axiom of our axiomatic approach. To 
this end we have to determine how individual packets of charge change in time as they move 
with velocity u % through space. A convenient way to describe this change is provided by the 
material derivative D/Dt which also is often called convective derivative [20]. It allows to 
calculate the change of a physical quantity as it appears to an observer or a probe that follows 
this quantity. Then electric charge conservation can be expressed as 

DO 

isr = - <8) 

where the material derivative is taken with respect to the velocity field u\ It can be rewritten 
in the following way 1201 . 



Dt Dt 



v{t) 
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^-dV + (f pu % ddi 
V(t) ot JdV(t) 

% + diifm')) dV . (9) 
<v(t) \ at 

Here we used in the last line the Stokes theorem in the form of (170b . The volume V(t) that is 
integrated over depends in general on time since it moves together with the electric charge that 
it contains. By means of ©, ©, and © we obtain the axiom of electric charge conservation 
in the local form as continuity equation, 

9 tP + ^J J = 0. (10) 

Now we use the inhomogeneous Maxwell equation © in order to replace within the con- 
tinuity equation (flOb the charge density by the divergence of D % . This yields 

dJdtD^ + r) =o. (ii) 



Again we invoke the Poincare lemma, now in the form (l68l) . and write the sum d t D l + J 1 as 
the curl of the integrand of a line integral which we denote by Hi. This yields 



e ijk d~H k - d t D l = J 1 



(curl H -V = J). (12) 



Equation (fT2b constitutes the remaining inhomogeneous Maxwell equation, the Ampere- 
Maxwell law, which, in this way, is derived from the axiom of charge conservation. The 
fields D % and Hi are called electric excitation (historically: electric displacement) and mag- 
netic excitation (historically: magnetic field), respectively. From (0) and <fT2l) it follows that 
their Si-units are [D 1 ] = As/m 2 and [Hi] = A/m. 

Some remarks are appropriate now: We first note that we obtain the excitations D % and 
Hi from the Poincare lemma and charge conservation, respectively, without introducing the 
concept of force. This is in contrast to other approaches that rely on the Coulomb and the 
Lorentz force laws Q. Furthermore, since electric charge conservation is valid not only on 
macroscopic scales but also in micropysics, the inhomogeneous Maxwell equations © and 
(fT2l are microphysical equations as long as the source terms p and J 1 are microscopically 
formulated as well. The same is valid for the excitations D % and H^ They are microphysical 
quantities — in contrast to what is often stated in textbooks, see 0, for example. We finally 
remark that the inhomogeneous Maxwell equations © and (fT2l can be straightforwardly put 
into a relativistically invariant form. This is not self-evident but suggested by electric charge 
conservation in the form of the continuity equation (TTOb since this fundamental equation can 
also be shown to be relativistically invariant. 



2.2 Lorentz force (axiom 2) and merging of electric and magnetic field 
strengths 

During the discovery of the electromagnetic field, the concept of force has played a major role. 
Electric and magnetic forces are directly accessible to experimental observation. Experimental 
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evidence shows that, in general, an electric charge is subject to a force if an electromagnetic 



field acts on it. For a point charge q at position x^ 
velocity u\ we postulate the Lorentz force 



we have p(x l ) = q5(x t — x q l ). If it has the 



Fi = q(Ei + e ijk u j B h 



(13) 



as second axiom. It introduces the electric field strength Ei and the magnetic field strength 
B l . The Lorentz force already yields a prescription of how to measure Ei and B % by means 
of the force that is experienced by an infinitesimally small test charge q which is either at 
rest or moving with velocity u\ Turning to the dimensions, we introduce voltage as "work per 
charge". In SI, it is measured in volt (V). Then [Fj]=VC/m and, according to (fT3t . [E l ] = V/m 
and [Bi] = Vs/m 2 = Wb/m 2 = T, with Wb as abbreviation for weber and T for tesla. 



q 







'B 1 1 




u 







Ei 





a) charge observed from its rest frame 



b) charge observed from inertial frame 
moving with respect to q 



Figure 1: A charge that is, in some inertial frame, at rest and is immersed in a purely magnetic 
field experiences no Lorentz force, see Fig. la. The fact that there is no Lorentz force should 
be independent of the choice of the inertial system that is used to observe the charge. There- 
fore, a compensating electric field accompanies the magnetic field if viewed from an inertial 
laboratory system which is in relative motion to the charge, see Fig. lb. 



From the axiom of the Lorentz force (fT31) . we can draw the conclusion that the electric and 
the magnetic field strengths are not independent of each other. The corresponding argument 
is based on the special relativity principle: According to the special relativity principle, the 
laws of physics are independent of the choice of an inertial system @. Different inertial 
systems move with constant velocities v l relative to each other. The outcome of a physical 
experiment, as expressed by an empirical law, has to be independent of the inertial system 
where the experiment takes place. 

Let us suppose a point charge q with a certain mass moves with velocity u % in an elec- 
tromagnetic field Ei and B\ The velocity and the electromagnetic field are measured in an 
inertial laboratory frame. The point charge can also be observed from its instantaneous iner- 
tial rest frame. If we denote quantities that are measured with respect to this rest frame by a 
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prime, i.e., by u'\ E'i, and B'\ then we have u n = 0. In the absence of an electric field in the 
laboratory system, i.e., if additionally E'i = 0, the charge experiences no Lorentz force and 
therefore no acceleration, 

F'i = q(E'i + e ijk u ,j B' k ) = . (14) 




Figure 2: The tetrahedron of the electromagnetic field. The electric and the magnetic excita- 
tions D\ Hi and the electric and the magnetic field strengths Ei, B l build up 4-dimensional 
quantities in spacetime. These four fields describe the electromagnetic field completely. Of 
electric nature are D l and Ei, of magnetic nature Hi and B\ The electric and the magnetic 
excitations D l , Hi are extensities, also called quantities (how much?), the electric and the 
magnetic field strengths Ei, B l are intensities, also called forces (how strong?). 



The fact that the charge experiences no acceleration is also true in the laboratory frame. 
This is a consequence of the special relativity principle or, more precisely, of the fact that the 
square of the acceleration can be shown to form a relativistic invariant. Consequently, 

F i = q{E i + e ijk vPB k ) = 0. (15) 

Thus, in the laboratory frame, electric and magnetic field are related by 

E l = -e l]k u^B k . (16) 

This situation is depicted in Fig[T] Accordingly, we find that electric and magnetic field 
strength cannot be viewed as independent quantities. They are connected to each other by 
transformations between different inertial systems. 

Let us pause for a moment and summarize: So far we have introduced the four electromag- 
netic field quantities D\ Hi and Ei, B l . These four quantities are interrelated by physical and 
mathematical properties. This is illustrated in Fig|2]by the "tetrahedron of the electromagnetic 
field". 
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2.3 Magnetic flux conservation (axiom 3) and the homogeneous Maxwell 
equations 

We digress for a moment and turn to hydrodynamics. Helmholtz was one of the first who 
studied rotational or "vortex" motion in hydrodynamics, see 0. He derived theorems for 
vortex lines. An important consequence of his work was the conclusion that vortex lines are 
conserved. They may move or change orientation but they are never spontaneously created nor 
annihilated. The vortex lines that pierce through a two-dimensional surface can be integrated 
over and yield a scalar quantity that is called circulation. The circulation in a perfect fluid, 
which satisfies certain conditions, is constant provided the loop enclosing the surface moves 
with the fluid 0. 

There are certainly fundamental differences between electromagnetism and hydrodynam- 
ics. But some suggestive analogies exist. A vortex line in hydrodynamics seems analogous to 
a magnetic flux line. The magnetic flux $ is determined from magnetic flux lines, represented 
by the magnetic field strength B\ that pierce through a two-dimensional surface S, 

$ := I B i ddi . (17) 
Js 

As the circulation in a perfect fluid is conserved, we can guess that, in a similar way, the 
magnetic flux may be conserved. Of course, the consequences of such an axiom have to be 
borne out by experiment. 

At first sight, one may find vortex lines of a fluid easier to visualize than magnetic flux 
lines. However, on a microscopic level, magnetic flux can occur in quanta. The corre- 
sponding magnetic flux unit is called flux quantum or fluxon and it carries $ — h/(2e) ~ 
2, 07 ■ 1CT 15 Wb, with h as Planck constant and e as elementary charge. Single quantized 
magnetic flux lines have been observed in the interior of type II superconductors if exposed 
to a sufficiently strong magnetic field, see [4J, p. 131. They even can be counted. The corre- 
sponding experiments provide good evidence that magnetic flux is a conserved quantity. 

But how can we formulate magnetic flux conservation mathematically? It is at this point 
instructive to reconsider the notion of the electric charge 

Q= [ pdv (18) 
Jv 

together with its corresponding conservation law 

d t Q+ [ J l da, t = 0. (19) 

JdV 

The rate of change of the electric charge within a specified volume V is balanced by the out- 
or inflowing charge across the surface dV. This charge transport is described by the electric 
charge current J % that is integrated over the enveloping surface dV. By means of the Stokes 
theorem in the form (l70l . equation (fT9l yields the local continuity equation 

d tP + d t J* = 0. (20) 
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Let us follow the same pattern in formulating magnetic flux conservation: Starting with the 
definition (ITvb of the magnetic flux, the corresponding conservation law, in analogy to ( fl9l) . 
reads 

<9 t $ + / J?dc* = 0, (21) 

JdS 

where we introduced the magnetic flux current Jf . Geometrically, this is a covariant vector 
that is integrated along a line OS, that is, along the curve bordering the 2-dimensional surface 
S. The conservation law (I2T1) tells us that the rate of change of the magnetic flux within 
a specified area S is balanced by the magnetic flux current Jf that is integrated along the 
boundary dS. Then the Stokes theorem in the form (1711) yields the local continuity equation 

+ e ijk djjf = 0. (22) 



One interesting consequence is the following: The divergence of d22l reads 

d l (d t B l ) = d i B i = p mae , <9 t p mag = 0. (23) 

Thus, we find a time-/«dependent term p mag , which acquires tentatively the meaning of a 
magnetic charge density. Let us choose a specific reference system in which p mag is constant in 
time, i.e., <9 t p mag = 0. Now we go over to an arbitrary reference system with time coordinate 
t 1 and spatial coordinates x % . Clearly, in general dt>p mag ^ 0. The only way to evade a 
contradiction to d23l is to require p mag = 0, that is, the magnetic field strength B % has no 
sources, its divergence vanishes: 



d i B i = 



(divB = 0). (24) 



This is recognized as one of the homogeneous Maxwell equations. Note that our derivation of 
(|24l was done under the assumption of magnetic flux conservation d2"T1) . Under this condition 
we find p mag = 0. 

In order to understand better the magnetic flux current, we note that Jf, as a covariant 
vector, has the same geometric properties as the electric field strength Ei. Additionally, Jf 
and Ei share the same physical dimension voltage/length, i.e., in SI, V/m. Accordingly, it is 
plausible to identify both quantities, 

jf = E t . (25) 

That also the sign chosen is the appropriate one (consistent with the Lenz rule) was discussed 
in 0. Then the local continuity equation d22l assumes the form 



d t B l + e ij %E k = 



(B + curl E = 0) . (26) 



This equation reflects magnetic flux conservation, the third axiom of our axiomatic approach. 
It also constitutes the remaining homogeneous Maxwell equation, that is, Faraday's induction 
law. 

At this point one might wonder to what extend the identification (l25l is mandatory. It turns 
out that it is special relativity that dictates this identification. We illustrate this circumstance 
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as follows: In the rest frame of a magnetic flux line B" the magnetic flux current vanishes, 
J'f = 0. The rest frame is also defined via the Lorentz force: In the absence of an electric 
field, E'i = 0, a test charge q is not accelerated by B n . Then a Lorentz transformation, 
together with (l24l) . yields an equation that relates B % and Jf in a laboratory frame according 
to 

jf = -t ijk u ] B k . (27) 

A comparison with (fToT) . which was obtained by an analogous transformation of a magnetic 
flux line from its rest frame to a laboratory frame, shows that the identification d25t needs to 
be valid, indeed. However, one should be aware that our simple argument requires E'i = in 
the rest frame of the considered magnetic flux line. 



2.4 Constitutive relations (axiom 4) and the properties of spacetime 

So far we have introduced 4x3 = 12 unknown electromagnetic field components D l ,Hi, 
Ei, and B l . These components have to fulfill the Maxwell equations ©, (fl"2l) . d24l . and d26l) . 
which represent 1 + 3 + 1 + 3 = 8 partial differential equations. In fact, among the Maxwell 
equations, only (fT2l and (l26b contain time derivatives and are dynamical. The remaining equa- 
tions, © and d24l . are so-called "constraints". They are, by virtue of the dynamical Maxwell 
equations, fulfilled at all times if fulfilled at one time. It follows that they don't contain in- 
formation on the time evolution of the electromagnetic field. Therefore, we arrive at only 
6 dynamical equations for 12 unknown field components. To make the Maxwell equations 
a determined set of partial differential equations, we still have to introduce additionally the 
so-called "constitutive relations" between the excitations D\ Hi and the field strengths E { , 
B\ 

The simplest case to begin with is to find constitutive relations for the case of electromag- 
netic fields in vacuum. There are guiding principles that limit their structure. We demand 
that constitutive relations in vacuum are invariant under translation and rotation, furthermore 
they should be local and linear, i.e., they should connect fields at the same position and at the 
same time. Finally, in vacuum the constitutive relations should not mix electric and magnetic 
properties. These features characterize the vacuum and not the electromagnetic field itself. 
We will not be able to prove them but postulate them as fourth axiom. 

If we want to relate the field strengths and the excitations we have to remind ourselves 
that E^ Hi are natural integrands of line integrals and D\ B l are natural integrands of surface 
integrals. Therefore, E if Hi transform under a change of coordinates as covariant vectors 
while D\ B l transform as contravariant vector densities. To compensate these differences 
we will have to introduce a symmetric metric field gij = gji. The metric tensor determines 
spatial distances and introduces the notion of orthogonality. The determinant of the metric 
is denoted by g. It follows that Jgg 1 ^ transforms like a density and maps a covariant vector 
into a contravariant vector density. We then take as fourth axiom the constitutive equations for 
vacuum, 

' ' (28) 



D* = e Qy /gcPE. 



3 1 
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H i =( f i 0> /g)- 1 g ij B>. 



(29) 



In flat spacetime and in cartesian coordinates, we have g = 1, g n = 1, and # u = for 
i ^ j. We recognize the familiar vaccum relations between field strengths and excitations. 
The electric constant e and the magnetic constant // characterize the vacuum. They acquire 
the Si-units [e ] = As/Vm and [/x ] = Vs/Am. 

What seems to be conceptually important about the constitutive equations (l28l) . d29l is that 
they not only provide relations between the excitations D\ Hi and the field strengths E^ B\ 
but also connect the electromagnetic field to the structure of spacetime, which here is repre- 
sented by the metric tensor g^. The formulation of the first three axioms that were presented 
in the previous sections does not require information on this metric structure. The connection 
between the electromagnetic field and spacetime, as expressed by the constitutive equations, 
indicates that physical fields and spacetime are not independent of each other. The constitutive 
equations might suggest the point of view that the structure of spacetime determines the struc- 
ture of the electromagnetic field. However, one should be aware that the opposite conclusion 
has a better truth value: It can be shown that the propagation properties of the electromagnetic 
field determine the metric structure of spacetime (4l|9|. 

Constitutive equations in matter usually assume a more complicated form than d28t . d29l) . 
In this case it would be appropriate to derive the constitutive equations, after an averaging 
procedure, from a microscopic model of matter. Such procedures are the subject of solid state 
or plasma physics, for example. A discussion of these subjects is out of the scope of this 
paper but, without going into details, we quote the constitutive relations of a general linear 
magnetoelectric medium: 



This formulation is due to Hehl & Obukhov |0]|5l[T5]|, an equivalent formulation of a "bian- 
isotropic medium" — this is the same as what we call general linear medium — was given 
by Lindell & Olyslager ifTTl [Toll . Both matrices e 1 ^ and /i^ 1 are symmetric and possess 6 
independent components each, is called permittivity tensor and /i^ 1 impermeability ten- 
sor (reciprocal permeability tensor). The magnetoelectric cross-term Yj, which is tracefree, 
7 fc fc = 0, has 8 independent components. It is related to the Fresnel-Fizeau effects. 

The 4-dimensional pseudo-scalar a, we call it axion piece J4|, represents one component. 
It corresponds to the perfect electromagnetic conductor (PEMC) of Lindell & Sihvola iTTA . a 
Tellegen type structure ll24ll25l . 

Accordingly, these pieces altogether, which we printed in d30b and (I3TT) in boldface for 
better visibility, add up to 6 + 6 + 8 + 1 = 20 + 1 = 21 independent components. The situation 
with 20 components is described in Post ifTBI (he reqiured a = without a real proof), that 
with 21 components in O'Dell ifToll . 

We can have 15 more components related to dissipation, which cannot be derived from a 
Lagrangian, the so-called skewon piece (see lfT4l and the literature given), namely 3 + 3 com- 
ponents of rik and m k (electric and magnetic Faraday effects), 8 components from the matrix 



Hi 



(e ij 



e llk n k ) Ej + ( + s/) B J + (a - a) B i , 
e ijk m k ) B j + (-7 J 'i + gj) Ej - (a + s) Ei . 



(30) 
(31) 
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(optical activity), which is traceless 5k k = 0, and 1 component from the 3-dimensional 
scalar s (spatially isotropic optical activity). This scalar was introduced by Nieves & Pal lfT3l . 
It has also been discussed in electromagnetic materials as chiral parameter, see Lindell et al. 
lfT2l . Note that s, in contrast to the 4-dimensional scalar a, is only a 3D scalar. We end then 
up with the general linear medium with 20 + 1 + 15 = 36 components. 

With the introduction of constitutive equations the axiomatic approach to classical elec- 
trodynamics is completed. We will see in the next Section [3]how this approach relates to the 
framework of gauge theory. 

3 On the relation between the axiomatics and the gauge ap- 
proach 

Modern descriptions of the fundamental interactions heavily rely on symmetry principles. In 
particular, this is true for the electromagnetic interaction which can be formulated as a gauge 
field theory that is based on a corresponding gauge symmetry. In a recent article this approach 
towards electromagnetism has been explained in some detail Q. The main steps were the 
following: 

• Accept the fact that physical matter fields (which represent electrons, for example) are 
described microscopically by complex wave functions. 

• Recognize that the absolute phase of these wave functions has no physical relevance. 
This arbitrariness of the absolute phase constitutes a one-dimensional rotational type 
symmetry U(i) (the circle group) that is the gauge symmetry of electromagnetism. 

• To derive observable physical quantities from the wave functions requires to define 
derivatives of wave functions in a way that is invariant under the gauge symmetry. The 
construction of such "gauge covariant" derivatives requires the introduction of gauge 
potentials. One gauge potential, the scalar potential 0, defines a gauge covariant deriva- 
tive Df with respect to time, while another gauge potential, the vector potential Ai, 
defines gauge covariant derivatives Df with respect to the three independent directions 
of space. 

• Finally, the gauge potentials </> and A, describe an electrodynamically non-trivial situa- 
tion, if their corresponding electric and magnetic field strengths 

Ei = -d l( p - d t Ai , (32) 
B i = e ijk djA k , (33) 

are non- vanishing. 

In the following we want to comment on the interrelation between the previously presented 
axiomatic approach and the gauge approach. It is interesting to see how the axioms find their 
proper place within the gauge approach. 



13 



3.1 Noether theorem and electric charge conservation 



In field theory there is a famous result which connects symmetries of laws of nature to con- 
served quantities. This is the Noether theorem which has been proven to be useful in both 
classical and quantum contexts. It is, in particular, discussed in books on classical electrody- 
namics, see fl9ll23ll . for example. 

Laws of nature, like in electrodynamics, e.g., can often (but not always) be characterized 
concisely by a Lagrangian density £ = £(^, d^, d t ^) which, in the standard case, is a 
function of the fields ^ of the theory and their first derivatives. Integration of the Lagrangian 
density £ over space yields the Lagrangian L, 



L = J £(^,d^,d t ^)dV, (34) 
and further integration over time yields the action S, 

S= I Ldt. (35) 



There are guiding principles that tell us how to obtain an appropriate Lagrangian density for 
a given theory. Once we have an appropriate Lagrangian density, we can derive conveniently 
the properties of the fields . For example, the equations of motion which determine the 
dynamics of \I> follow from extremization of the action S with respect to variations of ty, 

5^S = =>- equations of motion for . (36) 

Now we turn to the Noether theorem which connects the symmetry of a Lagrangian density 
£("$!, di^ff, d t ^) to conserved quantities. Suppose that £ is invariant under time translations <5 t - 
In daily life this assumption makes sense since we do not expect that the laws of nature change 
in time. Then the Noether theorem implies a local conservation law which expresses the con- 
servation of energy. Similarly, invariance under translations S x i in space implies conservation 
of momentum, while invariance under rotations 8 m j yields the conservation of angular mo- 
mentum, 

5 t £ = =>- conservation of energy , (37) 
S x i£ = =>- conservation of momentum , (38) 
5 u .j£ = =>- conservation of angular momentum . (39) 

These symmetries of spacetime are called external symmetries. But the Noether theorem also 
works for other types of symmetries, so-called internal ones — especially gauge symmetries. 
In this case, gauge invariance of the Lagrangian implies a conserved current with an associated 
charge. That is, if we denote a gauge transformation by 5 e we conclude 

5 e £ = =>- charge conservation . (40) 

If we apply this conclusion to electrodynamics, we have to specify the Lagrangian density to 
be the one of matter fields that represent electrically charged particles. Then invariance of 
this Lagrangian density under the gauge symmetry of electrodynamics yields the conservation 
of electric charge. Thus, if we accept the validity of the Lagrangian formalism, then we can 
arrive at electric charge conservation from gauge invariance via the Noether theorem. 
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3.2 Minimal coupling and the Lorentz force 



We already have mentioned that, according to (l36l) . we can derive the equations of motion of a 
physical theory from a Lagrangian density and its associated action. We can use this scheme to 
derive the equations of motion of electrically charged particles. In this case, the corresponding 
Lagrangian density (that of the electrically charged particles) has to be gauge invariant. 

If electrically charged particles are represented by their wave functions, the corresponding 
Lagrangian density will contain derivatives with respect to time and space. It follows that 
the Lagrangian density will be gauge invariant if we pass from partial derivatives to gauge 
covariant derivatives according to 

dt — Df:=d t + ^, (41) 
di — Df:=di-^Ai, (42) 

with q the electric charge of a particle, % = h/(2if) with h as the Planck constant and <p, Ai 
as electromagnetic potentials [3|. This enforcement of gauge invariance has a classical ana- 
logue. If electrically charged particles are represented by point particles, rather than by wave 
functions, we have to replace within the Lagrangian density the energy E and the momentum 
Pi of each particle according to [1231 

E — > E + q4>, (43) 
Pi — > Pi-qA. (44) 

The substitutions (|41|> . (|42|) or (l43l) . d4*4*h constitute the simplest way to ensure gauge invariance 
of the Lagrangian density of electrically charged particles. They constitute what commonly is 
called "minimal coupling". Due to minimal coupling, we relate electrically charged particles 
and the electromagnetic field in a natural way that is dictated by the requirement of gauge 
invariance. 

Having ensured gauge invariance of the action S, we can derive equations of motion by 
extremization, compare (l36b . It then turns out that these equations of motion contain the 
Lorentz force law (fT3l . Therefore the Lorentz force is a consequence of the minimal coupling 
procedure which couples electrically charged particles to the electromagnetic potentials and 
makes the Lagrangian gauge invariant. 

3.3 Bianchi identity and magnetic flux conservation 

The electromagnetic gauge potentials <fi and A4 are often introduced as mathematical tools to 
facilitate the integration of the Maxwell equations. Indeed, if we put the relations (l32l and (l33l 
into the homogeneous Maxwell equations (|24l and (I26T) . we recognize that the homogeneous 
Maxwell equations are fulfilled automatically. They become mere mathematical identities. 
This is an interesting observation since within the gauge approach the gauge potentials are 
fundamental physical quantities and are not only the outcome of a mathematical trick. Thus 
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we can state that the mathematical structure of the gauge potentials already implies the homo- 
geneous Maxwell equations and, in turn, magnetic flux conservation. In this light, magnetic 
flux conservation, within the gauge approach, appears as the consequence of a geometric iden- 
tity. This is in contrast to electric charge conservation that can be viewed as the consequence 
of gauge invariance, i.e., as the consequence of a physical symmetry. 

The mathematical identity that is reflected in the homogeneous Maxwell equations is a 
special case of a "Bianchi identity". Bianchi identities are the result of differentiating a po- 
tential twice. For example, in electrostatics the electric field strength E{ can be derived from 
a scalar potential according to 

Ei = d4 . (45) 
Differentiation reveals that the curl of Ei vanishes, 

& h d j E k = e iik d j d k <l> = Q, (46) 

which is due to the antisymmetry of e ljfe . Again, this equation is a mathematical identity, a 
simple example of a Bianchi identity. 

3.4 Gauge approach and constitutive relations 

The gauge approach towards electrodynamics deals with the properties of gauge fields, which 
represent the electromagnetic field, and with matter fields. It does not reflect properties of 
spacetime. In contrast to this, the constitutive equations do reflect properties of spacetime, 
as can be already seen from the constitutive equations of vacuum that involve the metric g^, 
compare (l28t and (l29t . Thus, also in the gauge approach the constitutive equations have to be 
postulated as an axiom in some way. One should note that, according to (I32L (l33t . the gauge 
potentials are directly related to the field strengths Ei and B\ The excitations D % and Hi are 
part of the inhomogeneous Maxwell equations which, within the gauge approach, are derived 
as equations of motion from an action principle, compare (l36b . Since the action itself involves 
the gauge potentials, one might wonder how it is possible to obtain equations of motion for the 
excitations rather than for the field strengths. The answer is that during the construction of the 
action from the gauge potentials the constitutive equations are already used, at least implicitly. 

4 Conclusion 

We have presented an axiomatic approach to classical electrodynamics in which the Maxwell 
equations are derived from the conservation of electric charge and magnetic flux. In the context 
of the derivation of the inhomogenous Maxwell equations, one introduces the electric and the 
magnetic excitation D % and Hi, respectively. The explicit calculation is rather simple because 
the continuity equation for electric charge is already relativistically invariant such that for the 
derivation of the inhomogeneous Maxwell equations no additional ingredients from special 
relativity are necessary. The situation is slightly more complicated for the derivation of the 
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homogeneous Maxwell equations from magnetic flux conservation since it is not immediately 
clear of how to formulate magnetic flux conservation in a relativistic invariant way. It should 
be mentioned that if the complete framework of relativity were available, the derivation of the 
axiomatic approach could be done with considerable more ease and elegance EJ. 

Finally, we would like to comment on a question that sometimes leads to controversial 
discussions, as summarized in 1 20], for example. This is the question of how the quantities 
Ei, D\ B\ and Hi should be grouped in pairs, i.e., the question of "which quantities belong 
together?". Some people like to form the pairs (Ei, B l ), (D\ Hi), while others prefer to build 
(Ei, Hi) , (D l , B % ). Already from a dimensional point of view, the answer to this question 
is obvious. Both, Ei and B % are vo/tage-related quantities, that is, related to the notions of 
force and work: In SI, we have [Ei] = V/m, [B l ] = T=Vs/m 2 , or [B l ] = [E { ] /velocity. 
Consequently, they belong together. Analogously, D l and Hi are current-related quantities: 
[D l ] = C/m 2 = As/m 2 , [Hi] = A/m, or [D % ] = [Hi] /velocity. Thermodynamically speaking, 
(Ei, B l ) are intensities (answer to the question: how strong?) and (D\ Hi) extensities (how 
much?) 

These conclusions are made irrefutible by relativity theory. Classical electrodynamics is a 
relativistic invariant theory and the implications of relativity have been proven to be correct on 
macro- and microscopic scales over and over again. And relativity tells us that the electromag- 
netic field strengths Ei, B % are inseparably intertwined by relativistic transformations, and the 
same is true for the electromagnetic excitations D\ Hi. In the spacetime of relativity theory, 
the pair (Ei, B l ) forms one single quantity, the tensor of electromagnetic field strength, while 
the pair (D l , Hi) forms another single quantity, the tensor of electromagnetic excitations. If 
compared to these facts, arguments in favor of the pairs (E iy Hi), namely that both are cov- 
ectors, and (D\ B % ), both are vector densities (see the tetrahedron in Fig. 2), turn out to be of 
secondary nature. Accordingly, there is no danger that the couples (Ei, B l ) and (D l , Hi) ever 
get divorced. 
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A Mathematical Background 

Within a theoretical formulation physical quantities are modeled as mathematical objects. The 
understanding and application of appropriate mathematics yields, in turn, the properties of 
physical quantities. In the development of the axiomatic approach, we made repeated use of 
integration, of the Poincare lemma, and of the Stokes theorem. It is with these mathemati- 
cal concepts that it is straightforward to derive the basics of electromagnetism from a small 
number of axioms. 



17 



A.l Integration 



Integration is an operation that yields coordinate independent values. It requires an inte- 
gration measure, the dimension of which depends on the type of region that is integrated 
over. We want to integrate over one-dimensional curves, two-dimensional surfaces, or three- 
dimensional volumes that are embedded in three-dimensional space. Therefore, we have to 
define line-, surface-, and volume-elements as integration measures. Then we can think of suit- 
able objects as integrands that can be integrated over to yield coordinate independent physical 
quantities. 



A.l.l Integration over a curve and covariant vectors as line integrands 

We consider a one-dimensional curve c = c(t) in three-dimensional space. In a specific 
coordinate system x\ with indices i = 1, 2, 3, a parametrization of c is given by the vector 

c(t) = (c\t),c 2 (t),c 3 (t)). (47) 

The functions c l (t) define the shape of the curve. For small changes of the parameter t, with 
t — > t + At, the difference vector between c(t + At) and c(t) is given by 

. . , /Ac 1 Ac 2 Ac 3 \ 

M*) = 77.77.77 (48) 



x At ' At ' At f 

compare Fig|3] In the limit where At becomes infinitesimally we obtain the line element 

dc(t) = {de l {t),d(?{t),dc*{t)) 

'dc^t) <9c 2 (t) dc\t)\ du (49) 



at ' dt ' at J 

It is characterized by an infinitesimal length and an orientation. 



We now construct objects that we can integrate over the curve c in order to obtain a co- 
ordinate invariant scalar. The line element dc contains three independent components dc l . If 
we shift from old coordinates x l to new coordinates yi' = y J ' (x l ) these components transform 
according to 

dc>' = Q -^dc l . (50) 

ox 1 

Therefore we can form an invariant expression if we introduce objects en = a(x l ), with three 
independent components a { , that transform in the opposite way, 

dx % 

w 



a 3' = 777 ft i • (51) 



This transformation behavior characterizes a vector or, more precisely, a covariant vector (a 
1-form). It follows that the expression 

«i dd = ay dc? (52) 



18 




yields the same value in each coordinate system. 

Thus, we can now immediately define integration over a curve by the expression 

J a.i dd = j ' ol\ dc 1 + «2 dc 2 + dc 3 
f ( dc 1 dc 2 dc 3 

The last line shows how to carry out explicitly the integration since a* and c l are functions of 
the parameter t. 




A.1.2 Integration over a surface and contravariant vector densities as surface inte- 
grands 

Now we consider a two-dimensional surface a = a(t, s). Within a specific coordinate sys- 
tem x l , a parametrization of a is of the form 

a(t, s) = (o 1 ^, s), a 2 (t, s), a 3 (t, s)) (54) 

with parameters t, s and functions a l (t, s) that define the shape of the surface. 



An elementary surface element is bound by lines t =const, t + dt =const, s =const, and 
s + ds =const, compare Fig El It is characterized by the two edges ^ dt and ^ ds. These 
edges span an infinitesimal surface, the area and orientation of which is characterized by a 
covariant vector da ri that points normal to the infinitesimal surface. The vector ddi is given by 
the vector product of ^ dt and ^ ds, 

. da? da k 

ddi = ejo-fe— — — dt ds . (55) 
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Figure 4: Parametrization of a surface a(t, s). The lines t =const, t + At =const, s =const, 
and s + As =const circumscribe a surface Acij that is spanned by the edges dt and ^ ds. 
In the limit At — > dt, As — > ds, it becomes an elementary surface element da,. 



In order to know how the components dai transform under coordinate transformations 
yj = yi we have to know the transformation behavior of the symbol e^. Since in any 
coordinate system, assumes the values 0, 1, or -1 by definition, it is obvious that in general 



dx 1 dx^ dx 
dy 1 ' dyi' dy 

This is because the determinant of the transformation matrix, i.e., 



doc^ doc^ dcc^ 

det{dx/dy)=e ijk ^-^-^, (57) 

is, in general, not equal to one. But it follows from (1571) that the correct transformation rule 
for 6ijk is given by 



e i'j'k' j_j./q_/q.a a_.ii a_.ii a_,k' €i i k 



dx 1 dxi dx k 
det(dx/dy) dy 1 ' dyi' dy k 

. ,c. ia .dx 1 dx j dx k 
= det{dy/dx)-^- 7 -^- ? -^e ijk . (58) 

With (1551 this yields the transformation rule for the components dai, 

dx 1 

day = det(dy / dx) —-jdcii . (59) 
oy J 

Now we construct quantities that can be integrated over a surface. Since a surface element 
is determined from three independent components dcii we introduce an integrand with three 
independent components f} 1 that transform according to 

^ det(dyfdx) dx 1 ^ ^ ^ 
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Transformation rules that involve the determinant of the transformation matrix characterize so- 
called densities. Densities are sensitive towards changes of the scale of elementary volumes. 
In physics they represent additive quantities, also called extensities, that describe how much 
of a quantity is distributed within a volume or over the surface of a volume. This is in contrast 
to intensities. The covariant vectors that we introduced as natural line integrals are intensive 
quantities that represent the strength of a physical field. 

The transformation behavior d60t of the components (3 % characterizes a contravariant vector 
density. With this transformation behavior the surface integral 

ftddi = / Few — gs dt ds (61) 
yields a scalar value that is coordinate independent. 



A.1.3 Integration over a volume and scalar densities as volume integrands 

We finally consider integration over a three-dimensional volume v in three-dimensional space. 
Again we choose a specific coordinate system x l and specify a parametrization of v by 

v(t, s, r) = (y x (t, s, r),v 2 (t, s, r),v 3 (t, s, r)) , (62) 

with three parameters t, s, and r. 

An elementary volume element dv is characterized by three edges ^ dt, ^jjds, and ^dr. 
The volume, which is spanned by these edges, is given by the determinant 

'dv 1 , dv 1 , dv 1 , 
dv = det ( -z—at, -rr—ds, -^—dr 
at as or 



dv 1 dyi dv k 
dt ds Or 



e ijk~^r~ — —dtdsdr. (63) 



It is not coordinate invariant but transforms under coordinate transformations yi' = yi' (x l ) 
according to 

dv' = det(dy/dx)dv. (64) 

Since the volume element dv constitutes one independent component, a natural object to 
integrate over a volume has one independent component as well. We denote such an integrand 
by 7. It transforms according to 

^ det(dy/dx) ^ ' 
This transformation rule characterizes a scalar density and yields 

, f dv 1 dv j dv k , 
'j dv = / 7 ^ijk-^-^—-^— dt ds dr (66) 



dt ds dr 



as a coordinate independent value. 
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A.2 Poincare Lemma 



The axiomatic approach takes advantage of the Poincare lemma. The Poincare lemma states 
under which conditions a mathematical object can be expressed in terms of a derivative, i.e., 
in terms of a potential. 

We consider integrands c^, (3\ and 7 of line-, surface-, and volume integrals, respectively, 
and assume that they are defined in an open and simply connected region of three-dimensional 
space. Then the Poincare lemma yields the following conclusions: 

1 . If oii is curl free, it can be written as the gradient of a scalar function /, 



ijk d jak = on = dj. (67) 



e 

2. If (5 l is divergence free, it can be written as the curl of the integrand ctj of a line integral, 

8^ = p = ^%a k . (68) 

3. The integrand 7 of a volume integral can be written as the divergence of an integrand (3 1 
of a surface integral, 

7 is a volume integrand 7 = di(3 l . (69) 

While conclusions (l67t . (l68t are familiar from elementary vector calculus, this might not be 
the case for conclusion (l69l . However, (1691 is rather trivial since, in cartesian coordinates 
x, y, z, for a given volume integrand 7 = 7(2;, the vector ft 1 with components (3 X = 
f*'y(t,y,z)/3dt, (3 V = /* 7(2:, t, z)/3 dt, and (3 Z = /* 7(ar, y, t)/3 dt fulfills Of course, 
the vector /3* is not uniquely determined from 7 since any divergence free vector field can be 
added to [3 l without changing 7. We further note that 7, as a volume integrand, constitutes a 
scalar density. It can be integrated as above to yield the components of (3 l as components of 
a contravariant vector density. Therefore the integration does not yield a coordinate invariant 
scalar such that 7 cannot be considered as a natural integrand of a line integral. 

A.3 Stokes Theorem 

In our notation Stokes theorem, if applied to line integrands 014 or surface integrands (3 l , yields 
the identities: 

I d i f3 i dv= I ftddi, (70) 

Jv JdV 

e ljk dja k dai= / aide 1 , (71) 
Jas 

where dV denotes the two-dimensional boundary of a simply connected volume V and OS 
denotes the one-dimensional boundary of a simply connected surface S. 
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